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DEFLAGRATION  TO  DETONATION  IN  THERMONUCLEAR  SUPERNOVAE 


1.  Introduction 

The  standard  model  for  Type  la  supernovae  (SN  la)  is  an  explosion  by  a  process  of 
thermonuclear  combustion.  One-dimensional  calculations  show  that  the  two  possible  modes 
of  this  combustion  are  either  supersonic  detonation  or  subsonic  deflagration  (Arnett  1969; 
Hansen  k  Wheeler  1969;  Ivanova  et  al.  1974;  Nomoto  et  al.  1976;  Nomoto  et  al.  1984). 
There  is  evidence,  both  theoretical  and  empirical,  that  SN  la  involve  a  transition  from  a 
deflagration  to  a  detonation,  or  DDT  (Khokhlov  1991a,b;  Arnett  k  Livne  1994;  Hoflich, 
Khokhlov  k  Wheeler  1995;  Hoflich  1995;  Wheeler  et  al.  1995;  Hoflich  k  Khokhlov  1996 
and  references  therein).  Pure  detonation  models  tend  to  produce  an  excess  of  iron  peak 
elements,  rather  than  the  intermediate  mass  elements  that  are  observed  in  the  outer  layers 
of  SN  la.  Models  in  which  DDT  occurs  can  account  for  the  properties  of  SN  la,  including 
the  observed  correlated  dispersion  in  light  curve  amplitude  and  decay  rate  (Hoflich  et  al. 
1995;  Hoflich  k  Khokhlov  1996).  In  such  models,  the  burning  front  starts  as  a  deflagration 
that  makes  a  transition  to  a  detonation  directly  (the  standard  “delayed  detonation” 
model),  or  after  the  star  has  expanded  due  to  partial  burning  and  then  recompressed  (the 
“pulsation  delayed  detonation”).  In  these  models,  the  density  at  which  DDT  occurs  is  a 
free  parameter.  The  comparison  of  these  models  with  observations  suggests  that  DDT  must 
occur  at  relatively  low  densities,  of  the  order  of  10’’  g  cm"^  .  In  this  paper  we  address  some 
of  the  key  physical  issues  of  DDT.  In  particular,  we  address  the  physical  mechanism  by 
which  DDT  may  occur  in  supernovae,  and  why  it  might  occur  at  low  densities.  The  process 
of  DDT  is  not  well  understood  in  either  supernovae  or  in  terrestrial  combustion.  Many  of 
the  principles  presented  here  apply  to  terrestrial  conditions  as  well. 

The  analysis  we  present  is  based  on  a  theory  of  unconfined  DDT  (Khokhlov,  Oran 
k  Wheeler  1996b).  It  is  assumed  that  the  inherent  mechanism  of  DDT  is  the  Zeldovich 
gradient  mechanism  (Zeldovich  et  al.  1970)  in  which  ignition  in  a  region  of  induction  time 
Manuscript  approved  October  10,  1996 
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gradient  leads  to  the  development  of  a  reaction  zone  and  a  shock  that  coalesce  to  form  a 
detonation.  In  this  theory,  it  is  presumed  that  a  mixed  region  with  gradients  of  temperature 
and  composition,  a  combination  of  which  results  in  an  induction  time  gradient,  is  created 
within  a  deflagration  region  by  turbulence  that  microscopically  mixes  hot  products  and  cold 
fuel.  Two  quantities  must  be  determined:  the  size  of  the  critical  region,  Lc,  required  for  a 
self-propagating  detonation  to  develop,  and  the  amount  of  mixing  (the  intensity  or  velocity 
of  turbulence)  required  to  produce  that  critical  length  of  mixing. 

In  supernovae,  there  are  two  situations  in  which  DDT  could  occur.  First,  DDT  could 
occur  directly,  when  turbulent  deflagration  alone  creates  a  large  mixed  region  which  ignites. 
For  this  process  to  occur,  the  flame  surface  associated  with  laminar  burning  has  to  be 
destroyed  by  turbulence.  When  the  laminar  flame  front  is  destroyed,  reactions  will  still 
continue  at  rates  depending  on  the  density  and  temperature,  but  microscopic  chemical  and 
thermal  mixing  of  fuel  and  product  occurs.  The  second  way  to  initiate  DDT  is  for  the 
deflagration  to  die  due  to  expansion  before  the  DDT  occurs.  If  insufficient  energy  has  been 
liberated,  the  star  will  expand  and  then  contract.  Mixing  of  fuel  and  products,  and  the 
subsequent  DDT  may  occur  during  the  contraction  phase. 

In  this  paper,  we  first  examine  the  mixing  process  and  determine  the  thermodynamic 
properties  of  the  mixture  as  a  function  of  the  degree  of  mixing.  We  show  that  these 
thermodynamic  properties  can  be  characterized  by  a  single  parameter  (e.g.,  the  mass 
fraction  of  fuel)  and  that  the  history  of  mixing,  that  is,  whether  mixing  occurs  before,  after, 
or  during  a  pulsation  phase  does  not  change  the  thermodynamic  properties  of  the  mixture 
significantly.  We  derive  the  critical  length  scale  to  initiate  a  detonation  by  the  Zeldovich 
mechanism.  Next,  we  evaluate  the  turbulent  intensity  necessary  to  break  a  laminar  flame 
and  to  create  the  conditions  for  DDT.  This  required  the  calculation  of  the  laminar  flame 
speed  and  laminar  flame  thickness  over  a  wide  range  of  densities,  3  x  10^  -  10^°  g  cm"^  . 
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Finally,  we  estimate  the  critical  density  for  DDT  based  on  the  existing  three-dimensional 
simulations  of  turbulent  deflagration  in  Type  la  supernovae,  and  discuss  implications  of  the 
results  for  supernovae  models. 


2.  Thermodynamics  of  mixing 

Consider  a  subsonic  mixing  process  that  occurs  at  approximately  constant  pressure, 
and  assume  that  no  reactions  occur  during  mixing.  The  conservation  of  mass  and  enthalpy 
during  mixing  gives  the  thermodynamic  parameters  of  the  mixture, 

^(Tm,  Pm)  =  Xf  h{Tf,  pf)  -I-  (1  —  Xf)  h{Tp,  Pp)  ,  (1) 

1  Xf  {1-  Xf) 

—  =  —  +  - - -  ,  (2) 

Pm  Pf  Pp 

where  h  is  the  enthalpy;  subscripts  f,  p,  and  m  stand  for  fuel,  product,  and  mixture  of 
products  and  fuel,  respectively;  Xf  is  the  fuel  mass  fraction  in  the  mixture.  We  assume 
that  the  state  of  the  products,  p,  is  determined  by  conservation  of  enthalpy  during  constant 
pressure  burning  in  a  deflagration  wave. 


h{Tp,  Pp)  —  h(Tf,  Pf)  Ag  , 

(3) 

P{Tp,Pp)  =  P{Tf,  Pf)  =  const  , 

(4) 

where  Ag  is  the  energy  released  by  nuclear  reactions. 

The  value  of  Ag  depends  on  the  density  at  which  the  laminar  flame  is  either  quenched 
by  turbulence  or  quenched  by  expansion,  and  is  not  known  exactly.  If  the  final  stage  of 
deflagration  takes  place  at  very  low  densities,  say,  less  than  lO^g  cm“^,  only  carbon  would 
be  able  to  react  in  a  deflagration  wave.  Temperatures  resulting  from  carbon  burning  would 
be  so  low  that  oxygen  would  not  react.  In  this  case,  we  expect  Ag  3.7  x  10^^  ergs  g“^, 
and  the  products  of  burning  would  consist  predominantly  of  half  silicon  and  half  oxygen  by 
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mass.  At  higher  densities,  both  carbon  and  oxygen  will  burn  to  form  silicon  group  elements. 
For  the  latter  case,  we  pick  a  typical  value  of  Ag  ~  6  x  10^^  ergs  g“^  Below  we  refer  to 
these  two  choices  as  Case  A  and  B,  respectively. 

Now  consider  mixing  and  changes  in  pressure  due  to  global  adiabatic  expansion  or 
contraction.  For  a  perfect  gas,  the  order  of  mixing  and  expansion  or  contraction  does  not 
change  the  final  state  of  mixed  material.  For  the  equation  of  state  of  white  dwarf  material, 
this  is  not  an  exact  statement,  but  we  show  that  it  is  a  good  approximation. 

Figure  1  shows  the  results  of  the  following  numerical  experiment  using  the  degenerate 
matter  equation  of  state  described  in  the  next  section.  A  half  carbon,  half  oxygen  fuel  at 
density  5  x  10®  g  cm“®  (point  /  in  Figure  1)  was  burned  at  constant  pressure  to  a  half  silicon 
and  half  oxygen  mixture  (to  reach  point  p)  and  then  mixed  with  X/  =  0.5  at  constant 
pressure  to  yield  the  conditions  at  point  mo-  Alternatively,  both  fuel  and  products  were 
expanded  adiabatically  by  a  factor  of  100  (to  points  /',  p'),  and  then  mixed  at  constant 
pressure  to  produce  the  conditions  at  point  m'.  Finally,  the  mixture  was  returned  along  the 
adiabat  to  the  initial  pressure.  The  final  state,  point  m,  is  very  close  to  trio,  even  imder 
this  extreme  adiabatic  expansion  and  contraction.  The  difference  in  temperature  between 
points  m  and  m'  is  5%,  and  the  difference  in  density  is  7%.  We  conclude  that,  to  a  very 
good  approximation,  the  thermodynamic  properties  of  a  mixture  do  not  depend  on  the 
history  of  the  mixing  and  expansion  processes.  We  make  that  assumption  in  the  following 
analysis  and  neglect  the  order  in  which  mixing  and  expansion  or  contraction  occur. 

Mixing  cold  fuel  and  hot  products  raises  the  entropy  of  the  mixture,  so  that  the 
mixture  will  ignite  at  lower  density  than  the  fuel  alone.  For  cases  A  and  B,  Figure  2  gives 
the  density  at  which  the  energy-release  rate  in  the  fuel-product  mixture  is  one-tenth  of  the 
maximum  energy-release  rate  for  isochoric  burning  at  this  density.  The  maximum  release 
rate  occurs  after  ignition,  as  the  temperature  rises  and  before  fuel  is  consumed.  Once  the 
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mixture  is  at  this  density,  the  induction  time  is  very  short,  only  ten  times  the  characteristic 
time  required  to  burn  all  of  the  fuel.  The  figure  shows  that  for  a  pure  fuel  to  ignite,  it 
needs  to  be  compressed  to  very  high  densities.  However,  even  for  fuel  mixed  with  rather 
small  amounts  of  hot  products,  the  ignition  density  decreases  drastically,  and  is  10®  - 10^ 
g  cm~®  for  a  wide  range  in  fuel  fraction.  If  matter  is  returning  to  high  densities  after  a 
pulsation,  a  mixture  containing  a  noticeable  amount  of  products  cannot  be  compressed  to 
densities  significantly  in  excess  of  10^  g  cm~®  before  it  ignites.  For  the  pulsation  delayed 
detonation  models,  therefore,  the  thermodynamic  considerations  place  a  limit  on  the 
conditions  at  which  DDT  can  occur:  DDT  can  only  occur  at  densities  near  or  below  10^  g 
cm  .  For  the  standard  delayed  detonation  models,  however,  no  such  limit  can  be  placed  by 
thermodynamic  considerations.  For  both  cases,  the  kinetics  and  hydrodynamics  of  ignition 
need  to  be  considered  to  determine  whether  DDT  can  occur  at  thermodynamically  allowed 
conditions. 


3.  Critical  length  scale  for  DDT 

We  now  determine  the  critical  size  of  the  mixed  region  capable  of  triggering  a 
detonation.  It  has  been  shown  (Khokhlov  et  al.  1996b)  that  there  are  two  processes  that 
are  necessary  for  DDT  to  occur  by  the  Zeldovich  mechanism.  First,  a  spontaneous  wave 
must  be  generated  and  a  shock-reaction  complex  must  be  formed.  This  complex  forms 
somewhere  in  the  middle  of  the  mixed  region  where  there  is  both  enough  fuel  and  high 
enough  temperature.  Second,  this  complex  must  be  able  to  survive  the  propagation  into 
the  unburned  material. 

To  estimate  the  critical  size  of  this  region,  Lc,  we  consider,  as  in  (Khokhlov  et  al. 
1996b),  a  nonuniform  region  created  by  mixing  the  products  of  isobaric  burning  and  fresh 
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fuel,  such  that  there  is  a  linear  spatial  distribution  of  fuel  fraction  Xf, 

( xtL,  0  <  X  <  L 

Xj{x)  =  ,  (5) 

(1,  x>  L 

where  L  is  the  size  of  the  mixed  region.  This  creates  a  region  of  oppositely  directed 
gradients  in  temperature  and  concentration  of  reactants.  Initially,  the  velocity  of  the 
material  is  zero,  and  the  pressure  Pq  is  constant  everywhere.  The  boundary  conditions  at 
X  =  0  are  reflecting  walls  (symmetry  conditions). 

The  system  is  described  by  the  Euler  equations  coupled  with  the  nuclear  kinetic 
equations 


^  +  V.(pU)  =  0, 

(6) 

^  +  V(pUU)  +  VP  =  0, 

(7) 

^  +  V-((B+P)U)  =  0, 

(8) 

dy 

+  (U-V)3>  =  K, 

(9) 

where  the  total  energy  density  E  —  Et  +  pU^/2  +  pNaQ  •  3^  is  the  sum  of  the 
thermal,  kinetic,  and  nuclear  energies,  U  is  the  fluid  velocity,  Et  is  the  thermal  energy 
density,  y  is  the  vector  in  the  composition  space  of  mole  fractions  of  nuclei,  Q  is  the 
vector  of  corresponding  binding  energies  per  nucleon,  and  'Il{p,Et,y)  is  the  vector 
of  net  rates  of  change  of  nuclear  species  with  time;  Na  is  the  Avogadro  number. 

To  describe  thermonuclear  burning,  an  alpha-nuclei  reaction  network  is  used  with 
y  =  {"^He,  “^Ne,  S,  ^^Ar,  ‘*°Ca,  The  equation 

of  state  includes  contributions  from  Fermi  gases  of  electrons  and  positrons,  radiation,  and 
ions.  Further  details  of  the  input  physics  can  be  found  in  (Khokhlov  1993). 

The  system  of  equations  is  integrated  numerically  using  a  one-dimensional  version 
of  the  time-dependent,  compressible  fluid  code  ALLA  based  on  the  Piecewise  Parabolic 
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Method  (Collela  k  Woodward  1984;  Collela  k  Glaz  1985),  implemented  as  lagrangian 
step  plus  remap.  Details  of  the  implementation  are  given  in  (Khokhlov  1995;  Khokhlov 
et  al.  1996a).  The  nuclear  reactions  are  coupled  to  fluid  dynamics  by  time-step  splitting. 
The  kinetic  equations  are  integrated  together  with  the  equation  of  energy  conservation 
using  a  stiff  solver  with  adjustable  substeps  to  keep  the  accuracy  better  than  1%.  Most 
of  computations  were  done  using  a  grid  with  1024  equidistant  computational  cells.  The 
convergence  of  the  numerical  solutions  was  tested  for  selected  cases  by  varying  the  number 
of  cells  from  1024  to  4096.  Computations  were  carried  out  in  both  planar  and  spherical 
geometries. 

Computations  were  performed  assuming  that  the  products  of  burning  that  formed 
the  mixture  (state  p  in  the  thermodynamic  relations  given  in  equations  1  to  4)  are 
characterized  by  the  value  of  the  energy  released  by  the  deflagration  Ag  =  4  x  10^^  ergs 
g“^  or  A?  =  3  X  10^^  ergs  g~^  These  values  are  close  to  Case  A  of  §  2  and  have  been 
selected  to  test  the  sensitivity  of  the  results  to  Aq.  For  Ag  =  4  x  10^'^,  the  composition  of 
products  was  taken  to  be  pure  .  For  the  second  value,  the  composition  of  products  was 
taken  to  be  half  and  half  by  mass.  The  system  was  prepared  in  an  initial  state 
and  then  allowed  to  evolve  in  time,  first  until  ignition  takes  place,  then  until  the  formation 
(or  failure)  of  detonation,  and  finally  to  the  time  when  the  generated  detonation  or  shock 
leaves  the  computational  domain.  The  results  of  these  computations  are  not  sensitive  to 
the  assumed  energy  release  and  composition.  They  mainly  depend  on  the  density  at  which 
explosion  occurs,  and  on  the  size  of  the  explosion  region. 

Figure  3  is  an  example  of  a  computation  in  which  a  detonation  successfully  forms. 
The  initial  density  of  the  fuel  po  =  1-0  x  10^  g  cm~^  and  the  size  of  the  mixed  region  is 
L  =  5.0  X  10®  cm.  Figures  3a  through  3e  show  profiles  of  density,  temperature,  pressure, 
velocity,  and  carbon  mole  fraction,  respectively,  for  selected  times.  In  the  initial  state,  the 
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density  and  carbon  mole  fraction  of  the  mixture  increase,  and  the  temperature  decreases  as 
a  function  of  position  x.  Because  of  such  opposing  temperature  and  concentration  gradients, 
the  ignition  does  not  occur  first  at  the  origin,  x  =  0,  but  is  shifted  to  some  point  x  >  0 
where  the  energy  generation  rate  is  maximum.  Prom  this  point,  the  ignition  spreads  out 
supersonically  as  a  spontaneous  wave.  A  detailed  discussion  of  spontaneous  burning  is  given 
in  Khokhlov  et  al.  (1996b).  The  overpressure  and  material  velocity  in  the  spontaneous  wave 
grow  with  time  as  the  phase  velocity  of  the  burning  front  decreases.  When  the  speed  of  the 
spontaneous  wave  approaches  the  speed  of  sound,  a  shock  wave  emerges,  and  a  shock  wave 
-  reaction  complex  forms.  The  profiles  1  to  3  illustrate  the  spontaneous  wave  propagation. 
There  is  no  shock  present  at  this  stage.  Matter  is  continuously  compressed  and  accelerated 
in  the  process  of  burning,  and  expands  and  slows  down  afterwards.  The  shock  -  reaction 
complex  forms  between  the  profiles  3  and  4.  The  complex  continues  to  propagate  through 
the  matter  in  which  density  increases  and  temperature  decreases  (profiles  4  and  5),  but  the 
strength  of  the  complex  grows.  Finally,  the  complex  passes  into  the  fully  unburnt  cold  fuel 
in  the  form  of  a  slightly  overdriven  detonation  (profile  6),  and  begins  its  relaxation  to  a 
Chapman- Jouguet  state. 

Figures  3f  and  3g  show  the  nuclear  composition  for  the  times  corresponding  to  profiles 
2  and  7,  respectively.  These  figures  illustrate  that  at  densities  ^  lO^g  cm“^,  the  important 
reactions  are  -f  ,  ^He  ,  and  the  reactions  on  and  "^Mg  that 

transform  the  products  of  -1-  burning  into  the  -group  elements.  Reactions 
involving  oxygen  are  less  important.  The  initial  amount  of  oxygen  in  the  mixture  remains 
practically  intact.  The  energy  release  is  approximately  Aq  ~  (3.5  -  4)  x  10^^  ergs  g"^ 
consistent  with  the  value  of  the  energy  release  for  the  products  of  deflagration  adopted 
above,  and  with  the  assumptions  about  the  energy  release  made  in  §  2.  Figures  4a  and 
4b  show  a  successful  initiation  of  a  detonation  at  a  density  five  times  higher,  po  =  5  x  10^ 
g  cm“®  .  At  this  density,  burning  leads  to  higher  temperatures  due  to  a  decrease  in  the 
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specific  heat  of  degenerate  matter.  As  a  result,  oxygen  burns  to  silicon  rather  quickly 
(Fig.4b),  and  the  energy  of  the  explosion,  Aq  ~  7  x  10^7  ergs  g“^,  is  almost  two  times  larger 
than  in  the  previous  case  po  =  lO^g  cm~^. 

The  shock-reaction  complex  can  survive  propagation  down  the  temperature  gradient 
and  can  grow  into  a  fully  developed  detonation  only  if  conditions  in  front  of  the  complex 
change  slowly  enough,  that  is,  if  L  is  large  enough.  If  the  variations  of  preshock  density  and 
temperature  with  distance  are  too  fast,  the  reaction  zone  structure  is  unable  to  adjust  to 
the  decreasing  reaction  rates  behind  the  shock,  and  is  unable  to  compensate  for  the  energy 
required  to  shock  denser  fuel. 

Figures  5a  to  5f  illustrate  the  failure  of  the  initiation  of  a  detonation  at  po  =  lO^g  cm“^ 
when  the  size  of  the  mixed  region  was  20  percent  smaller,  L  =  4  x  10®  cm,  than  in  the 
case  presented  in  Fig.  3.  Here,  too,  the  spontaneous  wave,  and  the  shock  -  reaction 
complex  form  (profiles  1  to  3).  However,  the  complex  does  not  survive  the  propagation 
down  the  temperatiure  gradient.  The  shock  and  reaction  separate.  The  shock  passes  into 
the  unburned  matter  leaving  the  reaction  wave  behind.  It  is  the  condition  that  the  complex 
must  survive  the  propagation  down  the  temperature  gradient,  not  the  formation  of  the 
complex  alone,  that  determines  the  critical  size  of  the  mixed  region. 

Figure  6  summarizes  the  results  of  a  series  of  many  computations  performed  to 
determine  Lg,  the  critical  size  of  the  region  capable  of  triggering  a  detonation,  as  a  function 
of  fuel  density.  The  value  of  Lg  is  a  very  sensitive  function  of  density.  It  is  virtually 
independent  of  assumptions  made  about  the  energy  release  Aq  and  composition  of  the 
products  formed  during  the  preceeding  deflagration  stage.  The  critical  length  depends, 
however,  on  the  energy  release  during  the  explosion  itself.  As  was  mentioned  earlier,  the 
energy  release  increases  with  increasing  background  density,  which  makes  the  dependence  of 
Lc  on  density  very  strong.  The  values  of  Lg  obtained  in  spherical  and  planar  geometries  are 


-  10  - 


practically  the  same,  because  the  radius  of  curvature  of  a  spherical  front  is  approximately 
the  same  as  the  critical  length  which  is  much  larger  than  the  detonation  wave  thickness. 

The  simulations  show  that  the  spontaneous  wave  grows  to  an  appreciable  strength,  and 
a  shock  -  reaction  complex  forms  in  a  part  of  the  mixed  region  with  a  high  fuel  fraction, 
X/  ~  0.8  —  0.9.  Thus,  it  is  this  mixture  that  is  important  for  DDT.  As  shown  in  Figure  2, 
explosion  of  such  a  mixture  can  take  place  at  densities  <  lO^g  cm”^  (Fig.  2).  On  the  other 
hand,  as  shown  in  Figure  6,  the  critical  length  for  DDT  is  a  sensitive  function  of  density; 
it  increases  steeply  at  low  densities  and  it  is  of  order  3  x  10^  —  4  x  10^  cm  at  densities  of 
5  X  10®  —  10^  g  cm~®.  This  length  is  a  factor  of  10^  — 10®  larger  than  the  size  of  a  detonation 
reaction  zone.  However,  even  at  these  low  densities,  this  length  is  small  compared  to  the 
radius  of  a  white  dwarf.  The  length  scale  Lc  is  small  enough  to  be  formed  by  turbulent 
mixing  in  a  relatively  small  portion  of  the  supernova  ejecta  for  p  ^  5  x  10®  g  cm“®. 


4.  Conditions  on  turbulent  intensity 

Here  we  discuss  the  conditions  necessary  to  produce  a  mixed  region  of  size  Lg  for  both 
the  direct  and  the  pulsation  modes  of  delayed  detonation.  First  consider  DDT  in  the  direct 
mode.  A  fundamental  assumption  is  that  turbulence  must  destroy  the  flame  front  so  that 
burned  and  unburned  materials  can  mix  microscopically.  To  break  the  flame  front,  the 
turbulence  velocity  t4,  at  a  scale  comparable  to  the  thickness  of  the  laminar  front,  must 
exceed  the  normal  speed  of  the  laminar  flame, 

Ub^K  Se,  (10) 

where  Se  is  the  laminar  flame  speed,  Se  is  the  laminar  flame  thickness,  and  A'  ~  1  is 
a  coefficient  that  describes  the  ability  of  the  laminar  flame  to  survive  in  a  turbulent 
environment  (Khokhlov  et  al.  1996b). 
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The  mechanism  for  breaking  the  flame  is  flame  stretching,  which  includes  the  eflPects  of 
strain  and  curvature.  These  effects  depends  on  the  properties  of  the  particular  flame  and 
on  the  spectrum  of  the  turbulence.  For  terrestrial  flames,  the  coefficient  K  significantly 
exceeds  unity  according  to  both  theory  and  experiments  (Poinsot,  Candel  k  Trouve  1996; 
Roberts  et  al.  1993).  For  the  thermonuclear  flame,  the  value  of  K  is  not  known  exactly. 
The  analysis  of  the  nonlinear  stabilization  of  the  flames  in  a  turbulent  field  suggests  K  =  8 
(Khokhlov  1995).  The  numerical  simulations  of  the  flame  presented  in  (Khokhlov  1995)  are 
consistent  with  a  large  value  of  K.  For  the  subsequent  analysis  in  this  paper,  we  use  the 
two  values  K  =  1  and  K  =  8. 

Figure  7  shows  the  laminar  flame  speed  Si  and  thickness  6i  for  a  0.5C  +  0.50  mixture 
computed  for  the  range  of  densities  3  x  10^  —  10^°  g  cm“^  ,  as  described  in  the  appendix. 
In  the  same  figure,  we  show  Si  according  to  Timmes  and  Woosley  (1992)  who  computed 
the  laminar  flame  speed  at  high  densities,  p  >  2  x  10®  g  cm"®  .  Our  results  are  in  good 
agreement  with  theirs  at  these  densities.  At  low  densities  the  flame  speed  decreases  with 
density  more  rapidly  than  at  high  densities,  and  reaches  a  very  low  value  of  5/  3  x  10^ 

cm/s  at  p  ~  3  X  10^  g  cm"®  .  The  computed  values  of  St  give  C4  assuming  K  =  1.  The 
curve  C/®  =  8Si  shows  the  required  turbulent  intensity  if  K  =  8. 

To  relate  t4  to  the  intensity  of  turbulent  motions  on  larger  scales,  additional 
assumptions  are  required  about  the  spectrum  of  turbulent  motions  inside  the  white  dwarf 
during  the  explosion.  The  most  favorable  assumption  for  DDT  is  that  the  turbulent 
spectrum  is  Kolmogorov.  By  assuming  that  a  Kolmogorov  spectrum  is  established,  we 
ignore  the  time  delay  required  to  establish  the  steady-state  cascade  of  turbulent  energy,  and 
ignore  freeze-out  of  turbulence  due  to  the  expansion  of  the  star.  Under  the  assumption  of  a 
Kolmogorov  spectrum,  the  speed  of  turbulent  motions  at  the  scale  of  the  laminar  front  6i  is 
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related  to  the  speed  of  turbulent  motions  on  larger  scales  L>  6i  as 


During  the  deflagration  phase  of  the  explosion,  large  portions  of  the  star  will  be 
occupied  by  the  very  convoluted  flame.  To  be  specific,  we  pick  a  typical  value  L  =  L/  ~  10® 
cm  for  the  maximum  size  of  the  turbulent  flame  region  (Khokhlov  1995).  The  turbulent 
region  may  be  somewhat  smaller  in  the  beginning  of  the  explosion,  and  somewhat  larger 
at  the  end,  but  this  value  is  of  the  correct  order  of  magnitude.  If  we  assume  that  the 
deflagration  is  strong  enough  to  unbind  the  star,  we  have  to  assume  that  Lj  is  large 
enough.  The  large-scale  turbulence  in  a  supernova  is  driven  by  buoyancy.  If  we  neglect 
the  effects  of  freeze-out  due  to  expansion  of  the  star,  we  expect  Ui^  be  of  order  of  the 
Rayleigh-Taylor  characteristic  velocity  at  this  scale  Ul^  —  Ur-t  —  O.b^JgLf  where  g  is 
the  effective  gravitational  acceleration.  In  a  star  in  hydrostatic  equilibrium,  or  close  to  it, 
the  value  of  -JgLf  does  not  exceed  the  speed  of  sound.  In  equilibrium,  P/R  pg,  where 
R  is  of  order  the  stellar  radius.  Prom  this  it  follows  that  ~  yJP/p  ~  and  hence 
^gLj  <  as  for  Lf  <  R.  For  our  estimates  we  take  the  value  oi  Ulj  =  Ur-t  =  O.Scs.  The 
values  of  the  turbulent  velocity  at  the  scale  of  the  flame  front  estimated  from  equation  (11) 
using  these  values  is  plotted  in  Figure  7.  The  intersection  of  the  curve  Us,  with  the  curves 
Si  and  C/®  gives  estimates  of  the  transition  density  ptr  <  (4  —  8)  x  10^  g  cm“®  below  which 
DDT  may  occur.  We  emphasize  that  this  estimate  of  ptr  is  made  neglecting  any  time  delay 
in  establishing  the  Kolmogorov  spectrum,  and  neglecting  freeze-out  of  turbulence  due  to 
expansion  of  the  star. 

In  order  to  take  the  effect  of  freeze-out  into  account,  three-dimensional  simulations  of 
the  explosion  of  the  entire  star  are  required.  To  date,  only  one  such  simulation  has  been 
performed  (Khokhlov  1995),  and  the  analysis  presented  in  that  paper  suggests  that  the 
expansion  freezes  scales  exceeding  approximately  L  ~  10®  —  10^  cm,  and  that  the  turbulent 
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burning  velocity  at  these  scales  do  not  exceed  C/x,  ~  10®  -  10’’  cm/s.  If  we  take  from  this 
range  the  most  favorable  values  for  DDT,  Ul  =  10’  cm/s  and  L  =  10®  cm,  equation  (11) 
gives  another  estimate  of  the  turbulent  velocity  at  the  scale  of  the  flame  front,  which  is 
also  shown  in  Figure  7.  The  intersection  of  the  curve  with  the  curves  C/®  and  Si  gives, 
we  believe,  a  more  realistic  estimate  of  the  transition  density  ptr  <  (2  -  5)  x  10’  g  cm"®  . 

There  are  other  effects  that  may  influence  the  transition  density.  Self-turbulization 
of  the  flame  on  scales  less  than  Lf  might  be  possible,  as  is  preconditioning  of  a  mixture 
by  shocks  generated  by  the  turbulence  itself.  The  latter  effects  might  somewhat  facilitate 
DDT,  but  not  to  a  great  extent  (Khokhlov  et  al.  1996b).  Self-turbulization  due  to  the 
Landau- Darrieus  instability  is  limited  by  nonlinear  effects.  Shock  preconditioning  requires 
a  very  strong  shock  to  be  effective.  We  believe  that  there  are  two  major  unknown  factors 
in  estimating  the  value  of  ptr^  the  spectrum  of  turbulence  in  the  exploding  star,  and  the 
extent  to  which  the  flame  can  survive  turbulent  stretch  on  small  scales  (coefficient  K). 

Now  consider  the  case  in  which  the  conditions  for  DDT  do  not  occur  in  the  deflagration 
phase,  and  the  deflagration  is  extinguished  by  the  expansion.  In  this  case,  burning  is 
quenched  by  expansion  and  the  products  and  fuel  mbc  freely.  There  is  no  laminar  flame 
front  and  there  are  no  nuclear  reactions.  This  situation  is  still  Rayleigh-Taylor  unstable. 
Previous  estimates  (Khokhlov  1991a;  Arnett  &  Livne  1994)  have  shown  that  the  mixed 
R-T  region  may  be  on  the  order  of  10®  -  10®  cm,  which  is  much  larger  than  Lc  at 
p  ^  5  X  10®  —  lO’g  cm  ®.  With  a  Kolmogorov  cascade,  the  eddy  turnover  time  scales  as 
A®/®,  where  A  is  the  size  of  the  eddy.  When  the  star  is  compressed,  the  turbulent  motions 
on  all  scales  will  be  enhanced  and  mixing  should  proceed  vigorously  (this  is  opposite  to 
the  effect  of  freeze-out  expected  during  the  expansion  of  the  star).  It  is  thus  reasonable  to 
expect  microscopic  mixing  to  be  complete  before  the  contraction  phase.’  This  suggests  that 


’We  have  invoked  microscopic  mixing  in  this  model.  There  is  considerable  evidence  that 
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the  pulsation  mechanism  can  lead  rather  naturally  to  conditions  where  DDT  is  unavoidable 
in  the  recompression  phase. 


5.  Conclusions 

We  applied  a  theory  of  DDT  in  unconfined  conditions  (Khokhlov  et  al.  1996a)  to  DDT 
in  Type  la  supernovae.  The  two  basic  assumptions  of  this  theory  are:  1)  the  gradient 
mechanism  is  the  inherent  mechanism  that  leads  to  DDT  in  unconfined  conditions,  and 
2)  the  mechanism  for  preparing  the  gradient  of  induction  time  is  turbulent  mixing  that 
requires  breaking  or  quenching  the  flame  front. 

Using  a  series  of  numerical  simulations,  we  determined  the  minimum  critical  size  of 
the  mixed  region,  Lc,  required  for  DDT  (§  3).  We  And  that  Lci  though  much  larger  than 
the  detonation  wave  thickness,  is  at  least  two  to  three  orders  of  magnitude  smaller  than 
a  white  dwarf  radius  for  p  >  5  x  10®  -  10^  g  cm"®  .  Thus,  a  mixed  region  of  critical  size 
required  for  triggering  a  detonation  can,  in  principle,  be  formed  by  turbulent  mixing  in  a 
small  portion  of  a  supernova. 

There  are  two  possible  modes  of  detonation  formation  in  a  Type  la  supernova:  a  direct 

mixing  occurs  in  supernova  ejecta  but  is  not  always  microscopic.  Several  lines  of  evidence 
strongly  suggest  that  the  mixing  in  SN  1987A  was  macroscopic,  but  not  microscopic.  We 
note,  however,  that  the  physical  conditions  here  and  in  a  collapse  environment  like  SN  1987A 
are  different.  The  instability  leading  to  mixing  in  SN  1987A  is  expected  to  be  a  Richtmeyer- 
Meshkov  instability  in  which  the  density/pressure  inversion  is  only  temporarily  induced  at 
composition/density  boundaries  due  to  shock  passage.  In  the  thermonuclear  explosions  that 
we  consider  here,  the  Rayleigh-Taylor  instability  is  intrinsic  and  long  lasting.  There  should 
be  ample  time  to  mix  to  small  length  scales. 
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mode  when  DDT  occurs  inside  an  active  deflagration  front,  and  a  pulsation  mode  which  is 
possible  if  DDT  does  not  occur  directly,  and  if  the  energy  released  during  the  deflagration 
phase  is  not  enough  to  unbind  the  star.  Then  the  star  will  experience  a  global  pulsation, 
and  mixing  will  take  place  between  the  cold  fuel  and  ashes  of  the  dead  deflagration.  The 
requirements  on  the  intensity  of  turbulent  mixing  needed  to  form  a  mixed  region  of  critical 
size  Lc  are  different  in  these  two  cases. 

For  a  direct  mode  of  DDT  (delayed  detonation  model),  the  intensity  of  turbulence  on 
small  scales  of  the  order  of  the  laminar  flame  thickness  must  be  high  enough  to  stretch 
and  break  the  surface  of  the  laminar  flame  front  (§  4,  equation  10).  We  computed  the 
laminar  flame  velocity  and  thickness  in  a  wide  range  of  densities  3  x  10^  —  10^®  g  cm“®  , 
and  estimated  the  critical  intensity  of  turbulence  at  a  scale  of  a  laminar  flame  thickness 
required  to  break  the  flame. 

We  also  tried  to  estimate  the  value  of  the  transition  density  ptr  at  which  DDT  should 
occur  during  a  Type  la  supernova  explosion.  Using  available  information  about  the 
intensity  of  turbulence  inside  the  exploding  star,  we  conclude  that  a  realistic  estimate  of  the 
transition  density  for  the  direct  delayed  detonation  explosion  is  ptr  (2  —  5)  x  lO^g  cm“®. 
However,  even  under  the  most  favorable  and  unrealistic  assumptions  (§4),  ptr  does  not 
exceed  10®  g  cm“®  . 

In  the  pulsation  mode  (pulsating  delayed  detonation  model),  microscopic  mixing 
throughout  large  areas  can  occur  because  there  is  no  flame  front  separating  burned 
and  unbmrned  matter,  and,  thus,  fuel  and  products  can  mix  freely.  We  considered  the 
thermodynamics  of  the  mixing  process  in  a  supernova  (§  2)  to  demonstrate  that  the  order 
of  mixing  and  global  expansion  or  contraction  of  the  supernova  is  not  important.  We  also 
found  that  even  a  small  fraction  (i^  1%)  of  high-entropy  products  of  burning  mixed  with 
cold,  low-entropy  fuel  increases  the  entropy  of  the  mixture  such  that  it  will  ignite  at  very 
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low  densities.  We  conclude  that  in  this  mode,  the  portions  of  the  mixture  of  size  comparable 
to  or  exceeding  Lc  would  unavoidably  detonate  around  ptr  —  10^  g  cm“®. 

Despite  various  details  of  DDT  that  require  future  consideration,  our  main  conclusion 
is  that  DDT  in  Type  la  supernovae  should  occur  in  the  density  range  5  x  10®  to  a  few 
times  10^.  At  higher  densities,  10®  —  10®  g  cm“®  ,  either  the  required  tiurbulent  velocity 
cannot  be  reached  (direct  mode  of  the  delayed  detonation)  or  the  material  inevitably 
ignites  before  such  high  densities  are  reached  (pulsating  mode  of  the  delayed  detonation). 
At  lower  densities,  the  size  of  a  critical  region  becomes  too  large.  This  conclusion  agrees 
with  the  results  of  global  modeling  of  Type  la  supernovae,  where  the  transition  density,  ptD 
at  which  DDT  occurs  is  a  free  parameter.  Global  modeling  suggests  that,  in  order  to  fit 
observed  light  curves  and  spectra  of  SNIa,  the  transition  density  should  be  ptr  —  10^  g  cm“® 
(Khokhlov  1991a;  Woosley  &  Weaver  1994;  Hoflich  et  al.  1995;  Nomoto  et  al.  1995;  Hoflich 
1995;  Hoflich  h  Khokhlov  1996). 
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A.  Properties  of  laminar  flames  in  degenerate  C-O  mixtures 

The  velocity  of  a  laminar  flame  Si  is  computed  as  the  eigenvalue  of  the  differential 
equations  describing  the  structure  of  a  steady-state  burning  wave  (Zeldovich  et  al.  1985). 
We  assume  that  the  wave  propagates  steadily  and  that  burning  occurs  at  constant  pressure 
since  the  flame  speed  is  much  less  than  the  sound  speed.  The  parameters  of  the  fiiel  and 
the  products  of  burning  are  designated  below  by  subscripts  0  and  1,  respectively.  The  mass 
flux  is  constant  through  a  steady  wave,  Jm  =  pu  =  PoSi  =  p\U\.  The  equation  of  energy 
balance  inside  the  wave  is 

,  dT  d  (  dT\ 

•""'‘'S  ~  di  j  +  K 

where  cp  is  the  specific  heat  at  constant  pressure,  x  is  the  thermal  conductivity,  Na  is  the 
Avogadro  number,  Q  is  the  vector  of  binding  energies  of  nuclei,  TZ  is  the  vector  of  net 
reaction  rates  (Section  3),  and  x  is  the  distance  across  the  wave.  The  terms  on  the  right 
side  of  (A.l)  describe  the  temperature  changes  due  to  the  heat  conduction  and  nuclear 
reactions.  Nuclear  mole  fractions  y  are  changing  inside  the  wave  due  to  nuclear  reactions, 

r  dy 

(A2) 

Ion  diffusion  is  ignored  in  (A.2)  since,  in  a  degenerate  gas,  the  mean  free  paths  of  ions  are 
much  smaller  than  the  mean  free  paths  of  electrons  and  photons.  The  density  p(r,  3^)  inside 
the  wave  is  related  to  T  and  y  via  the  constant  pressure  condition 

P(p,r,3^)  =  Po.  (A3) 

The  boundary  conditions  for  Eqs.(A.l)  to  (A.3)  are 

dT 

at  X  =  -00  :  T  =  To,  p  =  po,  ^  =  3^0,  T  =  (^-4) 


and 
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The  points  x  =  ±oo  are  singular  points  of  equations  (A.1)-(A.3).  The  integral  curve  which 
passes  through  these  singular  points  represents  the  solution  for  a  steady  flame  front.  We 
designate  the  value  of  Jm  corresponding  to  this  integral  curve  as  J*. 

For  low  temperatures,  T  <  Tg,  where  Tg  is  some  value  less  than  2  x  lO^K,  we  can 
neglect  the  energy  generation  term  in  (A.l)  in  comparison  with  the  term  describing  the 
heat  transport.  At  these  temperatures,  (A.l)  can  be  rewritten  as 

^  ^  {A.1) 

dT  X 

where  W  =  x%-  The  corresponding  initial  condition  for  (A.7)  is  W(ro)  =  0.  The 
integration  of  (A.7)  from  To  to  gives  a  certain  value  of  We  —  W{Te),  from  which  we  And 
the  corresponding  value  of  (^)j,  •  Prom  this  point,  we  integrate  the  full  equation  (A.l) 
with  the  energy  generation  term  included. 

The  value  J*  is  found  using  the  method  of  trial  and  error  (Zeldovich  et  al.  1985). 
The  integral  curves  for  Jm  ^  J*  diverge  from  the  J*  curve  when  x  — >  oo.  For  Jm  <  J* 
the  asymptotic  (x  — ♦  oo)  behavior  is  T  —*  — oo,  ^  — oo.  For  Jm  >  J*  the  asymptotic 

behavior  is  T  ^  oo,  ^  ^  oo.  We  use  the  difference  in  the  behavior  of  the  integral 
curves  for  discriminating  between  the  cases  Jm  <  J*  and  Jm  >  J*-  The  next  trial  value 
of  Jm  is  found  as  the  average  of  the  two  most  recent  trial  values  lying  on  different  sides 
of  J*.  It  requires  ~  20  iterations  to  obtain  the  value  of  J*  with  an  accuracy  of  10“^. 
Then  the  laminar  flame  speed  is  calculated  as  Si  =  J*/po-  The  values  of  6i  are  estimated 
approximately  as  the  thickness  of  a  zone  where  most  of  the  nuclear  energy  is  released  and 
the  CO  mixture  is  converted  into  Si-peak  elements.  Computations  are  done  for  the  initial 
mixture  of  equal  masses  of  and  with  a  small  amount  of  added  in  order  to  get 
a  neutron  excess  rj  =  0.002.  The  initial  temperature  is  To  =  10®K,  and  we  used  Tg  =  10® 
K.  Flame  parameters  are  not  sensitive  to  variations  of  To  and  Te  if  these  temperatures  are 
below  2  X  10®K. 
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In  computations  of  Si  and  (5/,  the  following  133  species  reaction  network  is  used;  n, 
p,  ^He,  20.22 jVe,  ^^Na, 

36-41  Ar^  39-43^^  '>0-46(70,  43-475^,  44-50 47-52y^  48-56(;^^  “-62 /Tg^  55-61(7o, 

56-64^^^  57-65(7^^  forwaxd  reaction  rates  are  taken  from  the  Thielemann’s 

reaction  rate  library  (Thielemann,  Arnould  &  Truran  1987).  The  binding  energies  and 
partition  functions  were  taken  according  to  Woosley  et  al.  (1987).  The  backward  reaction 
rates  are  calculated  according  to  the  principle  of  detailed  balance.  Screening  corrections  are 
applied  to  the  most  important  reactions,  C,  ^2(7 +i6  i6q  ^le  q  3  4^g  _^i2 

which  control  the  thermonuclear  energy  release  at  low  temperatures.  Even  for  these 
reactions  screening  factors  were  found  to  be  relatively  small  (E  ^  2  -  5)  for  conditions  where 
nuclear  reactions  significantly  contribute  to  the  right  side  of  Eq.(2).  Most  of  the  nuclear 
reactions  occur  at  high  temperatures  and  the  corresponding  screening  factors  E  ^  1.5  are 
small  in  comparison  with  the  expected  uncertainties  in  the  nuclear  reaction  rates  (~  2  -  3). 
Screening  corrections  are  taken  from  Yakovlev  and  Shalybkov  (1989). 

The  total  thermal  conductivity  is  calculated  as  the  sum  of  electron  and  photon 
conductivities 

X  =  Xe  +  X7>  (>l-8) 

with  Xi  =  ^acT^/Snp.  The  expression  for  the  photon  opacity  (Paczynski  1983) 

86 

approximates  the  low  density  (;^  6  lO'^g  cm“^)  results  of  Buchler  and  Yueh  (1976)  that  take 
the  degeneracy  of  electrons  and  relativistic  scattering  into  account.  Yg  is  the  electron  mole 
fraction.  At  higher  densities,  the  contribution  of  the  photon  opacity  to  the  total  opacity  is 
very  small.  Our  results  at  high  densities  are  in  a  good  agreement  with  those  of  Timmes  & 
Woosley  (1992)  who  used  different  approximations  for  the  photon  thermal  conductivity. 

The  electron  thermal  conductivity  can  be  expressed  in  terms  of  the  effective  electron 


(1  +  2.7  io‘V/r2)  V  V4.5  10*/ 


-1 


cm^ 

g 


(A.9) 
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collisional  frequency,  Ue,  as  (Yakovlev  and  Urpin  1980) 


Xe  =  4.09  10®  T 


X'’ 


'lOi®' 


ergs 


^/T+1^  \  Ve  j  cm  s  K 


(.4,10) 


where  x  —  p^/meC  ~  1.009  10  Pf  is  the  electron  Fermi  momentum,  and  Vf.  is  a 

sum  of  ion-electron  and  electron-electron  collisional  frequencies. 


=  ^ei  +  J^ee- 


(All) 


The  electron  -  ion  collisional  frequency  is  calculated  using  the  mean  ion  approximation 
following  Urpin  &  Yakovlev  (1980) 

_ Y  7'^K 

Uei  =  1.78  10i®\/r+^  (A12) 

■*  e 

where  Y  =  J2Yi  and  Z  =  Y^jY  are  the  mean  ion  mole  fraction  and  the  mean  atomic 
number,  respectively,  and  A  is  the  Coulomb  logarithm  which  can  be  expressed  for  gaseous 
and  liquid  ionic  states  as  the  sum  of  the  logarithm  obtained  in  the  Born  approximation  and 
the  non-Born  correction  (Yakovlev  and  Urpin  1980,  Yakovlev  1987) 


A  =  Ab  +  5A, 


(.4,13) 


with 


and 


(A.14) 

(.4,15) 


TT  2  1  +  1.30a 

'  ^  2“^  l+a2(0.71- 0.54^02)  ’ 

a  X  ^  Z 

“  n/T+F  ’  137/?  ■ 


(.4.16) 

(A.n) 


The  thermal  conductivity  due  to  electron-electron  collisions  of  relativistic  electrons  in 
a  degenerate  electron  gas  was  calculated  in  the  static  screening  approximation  by  Urpin 


-  21  - 


and  Yakovlev  (1980), 


where 


(>1.18) 

(^.19) 


and 

T,  =  ^  =  3.307  10=  K  (>1.20) 

is  the  electron  plasma  temperature.  Recently  Yakovlev  (unpublished)  recalculated  the 
integral  entering  the  general  expression  of  numerically  and  obtained  the  following 
analytical  expression  for  J{y) 


Ay) 


1 

3 


(A21) 


with  a  =  0.113  and  b  =  1.247.  The  largest  error  in  the  fit  is  ~  6%  at  y  ~  2. 
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Aq  =  6.8  10’^  ergs/g 
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Fuel  Mass  Fraction 

Fig.  2  -  Ignition  density  as  a  function  of  the  fuel  fraction.  (Solid  line)  -  Aq 
ergs  g~^;  (Dashed  line)  -  Aq  =  6.8  x  10^^  ergs  g“^ 
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(a) 

Fig.  3  -  Initiation  of  a  detonation  for  po  =  10^  g  cm"^  and  L  =  10®  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10"®  (1),  7.05  x  10"®  (2), 
1.40  X  10-^  (3),  1.63  X  10-^  (4),  1.90  x  lO''*  (5),  2.33  x  10"^  (6),  2.76  x  10"^  (7)  and 
3.27  X  lO"'*  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon 
mole  fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g) 
Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  7. 
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(b) 

Fig.  3  (Continued)  —  Initiation  of  a  detonation  for  Po  =  10’  g  cm"’  and  L  =  10^  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10"^  (1),  7.05  x  10“^  (2) 

1.40  X  10-''  (3),  1.63  X  lO--*  (4),  1.90  x  10"''  (5),  2.33  x  lO--*  (6),  2.76  x  lO--*  (7),  and 
3.27  X  lO--*  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g)  Distribution  of 
nuclear  species  at  the  time  corresponding  to  profile  7. 
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(c) 

Fig.  3  (Continued)  —  Initiation  of  a  detonation  for  pq  =  10’  g  cm"’  and  L  =  lO’  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10"’  (1),  7.05  x  10“’  (2) 

1.40  X  lO"'*  (3),  1.63  X  10"^  (4),  1.90  X  10"^  (5),  2.33  X  lO"'*  (6),  2.76  x  lO"-*  (7),  and 
3.27  X  lO"-*  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g)  Distribution  of 
nuclear  species  at  the  time  corresponding  to  profile  7. 
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(d) 

Fig.  3  (Continued)  —  Initiation  of  a  detonation  for  po  =  10’  g  cm"^  and  I  =  10^  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10'*  (1),  7.05  x  10”^  (2), 

1.40  X  10"''  (3),  1.63  X  10-“  (4),  1.90  x  10”“  (5),  2.33  x  10"“  (6),  2.76  x  10-“  (7),  and 
3.27  X  10-“  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g)  Distribution  of 
nuclear  species  at  the  time  corresponding  to  profile  7. 
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(e) 

Fig.  3  (Continued)  —  Initiation  of  a  detonation  for  Po  =  10’  g  cm“^  and  L  =  10^  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10”*  (1),  7.05  x  10'^  (2), 

1.40  X  lO-^*  (3),  1.63  X  lO--*  (4),  1.90  x  10-“  (5),  2.33  X  10"“  (6),  2.76  X  10““  (7),  and 
3.27  X  10"“  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g)  Distribution  of 
nuclear  species  at  the  time  corresponding  to  profile  7. 
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Fig.  3  (Continued)  —  Initiation  of  a  detonation  for  Po  =  10’  g  cm“^  and  L  =  10^  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10"’  (1),  7.05  x  10“’  (2), 

1.40  X  lO"''  (3),  1.63  X  10"^  (4),  1.90  x  10"^  (5),  2.33  x  lO"-*  (6),  2.76  x  lO""  (7),  and 
3.27  X  10"“*  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g)  Distribution  of 
nuclear  species  at  the  time  corresponding  to  profile  7. 
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(g) 

Fig.  3  (Continued)  —  Initiation  of  a  detonation  for  Po  =  10’  g  cm"^  and  L  =  10^  cm.  Profiles 
are  shown  for  the  initial  time  t  =  0,  and  for  times  t  =  2.86  x  10’*  (1),  7.05  x  10“^  (2), 

1.40  X  lO-"  (3),  1.63  X  10-"  (4),  1.90  X  10'^  (5),  2.33  x  lO-"  (6),  2.76  X  10-“  (7),  and 
3.27  X  10"“  (8)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fi-action.  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  2.  (g)  Distribution  of 
nuclear  species  at  the  time  corresponding  to  profile  7. 
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Fig.  4  (Continued)  —  Initiation  of  a  detonation  for  Po  =  5  x  10’  g  cm~^  and 
L  =  1.3  X  0^  cm.  Profiles  are  shown  for  the  initial  time  t  =  0,  and  for  the  moments  of 
time  t  =  1.74  x  10“®  (1),  7.89  x  10"*  (2),  1.38  X  10”’  (3),  1.90  x  10”’  (4),  2.41  x  10”’  (5), 
3.38  X  10”’  (6),  4.31  x  10”’  (7)  and  4.96  x  10”’  (8)  sec.  (a)  Pressure,  (b)  Distribution  of  nuclear 
species  at  the  time  corresponding  to  profile  7. 
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(a) 

Fig.  5  -  Failed  initiation  of  a  detonation  for  po  =  1  x  10'^  g  cm“^  and  L  =  8  x  10^  cm. 
Profiles  are  shown  for  the  initial  time  t  =  0,  and  for  the  moments  of  time  t  =  5.64  x  10"® 
(1),  1.12  X  10-4  (2),  1.68  X  10-4  (3),  2.22  x  10-4  (4),  2.77  x  10-4  (5)  3  3^  ^  iq-4  (g) 

sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole  fraction,  (f) 
Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  5. 
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(b) 

Fig.  5  (Continued)  —  Failed  initiation  of  a  detonation  for  Po  =  1  x  10’  g  cm"^  and 
L  =  8  X  10''  cm.  Profiles  are  shown  for  the  initial  time  t  =  0,  and  for  the  moments  of  times 
r  =  5.64  X  10-*  (1),  1.12  X  10-'*  (2),  1.68  x  10“^  (3),  2.22  x  10'^  (4),  2.77  x  10'^  (5),  and 
3.31  X  10“''  (6)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  5. 
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(c) 

Fig.  5  (Continued)  —  Failed  initiation  of  a  detonation  for  Po  =  1  x  10’  g  cm“^  and 
L  =  8  X  10''  cm.  Profiles  are  shown  for  the  initial  time  r  =  0,  and  for  the  moments  of  times 
t  =  5.64  X  10“''  (1),  1.12  X  10-“  (2),  1.68  x  10"“  (3),  2.22  x  10““  (4),  2.77  x  10““  (5),  and 
3.31  X  10"“  (6)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  5. 
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Fig.  5  (Continued)  —  Failed  initiation  of  a  detonation  for  po  =  1  x  10’  g  cm"^  and 
L  =  8  X  10''  cm.  Profiles  are  shown  for  the  initial  time  t  =  0,  and  for  the  moments  of  times 
t  =  5.64  X  10“^  (1),  1.12  X  lO"'*  (2),  1.68  x  10“''  (3),  2.22  x  10“''  (4),  2.77  x  10“''  (5),  and 
3.31  X  10“^  (6)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  5. 
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Fig.  5  (Continued)  —  Failed  initiation  of  a  detonation  for  Po  =  1  x  10“'  g  cm"^  and 
I  =  8  X  10'*  cm.  Profiles  are  shown  for  the  initial  time  t  =  0,  and  for  the  moments  of  times 
t  =  5.64  X  10-=  (1),  1.12  X  lO-'*  (2),  1.68  x  10"''  (3),  2.22  X  10"^  (4),  2.77  x  10'^  (5),  and 
3.31  X  10-“'  (6)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  5. 
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(f) 

Fig.  5  (Continued)  —  Failed  initiation  of  a  detonation  for  Pp  =  1  x  lO’  g  cm~^  and 
8  X  10  cm.  Profiles  are  shown  for  the  initial  time  r  =  0,  and  for  the  moments  of  times 
r  =  5.64  X  10-^  (1),  1.12  X  10-''  (2),  1.68  x  10"''  (3),  2.22  x  10“^  (4),  2.77  x  10'^  (5).  and 
3.31  X  10-'*  (6)  sec.  (a)  Pressure,  (b)  Temperature,  (c)  Density,  (d)  Velocity,  (e)  Carbon  mole 
fraction,  (f)  Distribution  of  nuclear  species  at  the  time  corresponding  to  profile  5. 
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Fig.  6  -  Critical  length  for  the  initiation  of  a  detonation  as  a  function  of  the  fuel 
density.  (Circles)  -  Aq  =  4  x  10^^  ergs  g~^  for  deflagration  products.  (Triangles)  - 
Aq  =  3  X  10^^  ergs  g“^  for  deflagration  products  (all  but  one  triangle  coincide  with  the 
circles). 
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Fig.  7  Velocitiss  and  flame  thickness  as  a  function  of  density.  Si  and  5/  •“  laminar 
flame  speed  and  thickness,  computed  in  the  appendix.  S'f^  ^  laminar  flame  speed  according 
to  Timmes  &  Woosley  (1992).  -  critical  turbulent  velocity  assuming  =  8.  — 

estimated  turbulent  velocity  at  the  scale  of  the  flame  front  during  the  supernova  explosion, 
assuming  no  turbulence  freeze-out.  -  same,  but  assuming  turbulence  freeze-out. 


